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Abstract 



o 

\-{ ' Various aspects of relativistic calculations of the Sunyaev-Zeldovich effect are ex- 

c/3 . plored and clarified. We first formally show that the main previous approaches to 

^ ' the calculation of the relativistically generalized thermal component of the effect 

are equivalent. Our detailed description of the full effect results in a somewhat im- 
^ ■ proved formulation. Analytic approximations to the exact calculation of the change 

of the photon occupation number in the scattering, An, are extended to powers of 
the gas temperature and cluster velocity that are higher than in similar published 
treatments. For the purely thermal and purely kinematic components, we obtain 
identical terms up to the highest common orders in temperature and cluster ve- 
locity, and to second order in the Thomson optical depth, as reported in previous 
treatments, but we get slightly different expressions for the terms that depend on 
both the gas temperature and cluster velocity. We also obtain an accurate expression 
for the crossover frequency. 
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1 Introduction 



Compton scattering of the cosmic microwave background (CMB) by electrons 
in the hot gas in clusters of galaxies (Sunyaev & Zeldovich 1972) imprints 
on the radiation a spectral signature - the Sunyaev-Zeldovich (S-Z) effect - 
that has long been recognized as a uniquely important feature which can be 
used as an indispensable cosmological probe. Measurements with single-dish 
radio telescopes and interferometric arrays have led so far to the detection 
and imaging of the effect in more than 60 clusters. These measurements yield 
directly the properties of the hot intracluster (IC) gas, and indirectly the total 
mass of the cluster, insight on the evolution of clusters, and the values of the 
basic cosmological parameters (for reviews, see Rephaeli 1995a, Birkinshaw 
1999, Carlstrom et al. 2002). 

Exact description of the effect is needed in order to extract precise information 
from current measurements, and especially from upcoming high-frequency ob- 
servations. The calculations of Sunyaev & Zeldovich (1972) were based on a 
solution to the Kompaneets (1957) equation, a nonrelativistic diffusion ap- 
proximation to the exact kinetic (Boltzmann) equation that describes the 
scattering. The result of their treatment is a simple expression for the inten- 
sity change resulting from scattering of the CMB by electrons with a thermal 
velocity distribution. The effect has a second kinematic (Doppler) component 
which is proportional to the component of the cluster peculiar velocity along 
the hue of sight, Vr (Sunyaev & Zeldovich 1980). 

Rephaeli (1995b) has shown that the approximate nonrelativistic treatment 
of Sunyaev & Zeldovich (1972) is not sufficiently accurate for use of the effect 
as a precise cosmological probe. High electron velocities in the IC gas, and 
relatively large photon energy change in the scattering, necessitate a more 
exact relativistic calculation. Using the exact probability distribution, and the 
relativistically correct form of the electron Maxwellian velocity distribution, 
Rephaeli calculated the resulting intensity change in the limit of small op- 
tical depth to Thomson scattering, r, keeping terms linear in r. This semi- 
analytic calculation demonstrated that the relativistic spectral distribution 
of the intensity change is appreciably different from that derived by Sunyaev 
& Zeldovich (1972). Deviations from the Sunyaev & Zeldovich expression arc 
especially large near the crossover frequency, where the purely thermal ef- 
fect vanishes, and on the Wien side where the effect of boosting photons to 
relatively high frequencies by energetic electrons is important. 

The results of the calculations of Rephaeli (1995b) generated considerable 
interest that led to various generalizations and extensions of the relativistic 
treatment. Challinor & Lasenby (1998) obtained an analytic approximation to 
the solution of the relativistically generalized Kompaneets equation. Nozawa 
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et al. (1998a) improved the accuracy of this approximation by expanding to 
fifth order in = kTe/rrieC^, where Tg is the electron temperature. Sazonov 
& Sunyaev (1998) and Nozawa et ah (1998b) have extended the relativistic 
treatment also to the kinematic component obtaining - for the first time - 
the leading cross terms in the expression for the combined (thermal and kine- 
matic) intensity change, Alt + A/fe, which depends on both Tg and Vc, the 
cluster peculiar velocity. An analytic fit to the numerical solution, valid for 
0.02 < Ge < 0.05, and x = hu/kT < 20 (for t/ < 1130 GHz; T is the CMB 
temperature), has been given by Nozawa et al. (2000). 

Cluster X-ray and S-Z measurements have significantly improved in the last 
few years. Many nearby and distant clusters have been sensitively mapped 
spectrally and spatially by the Chandra and XMM satellites, yielding a wealth 
of information on the complex morphology and thermal structure of IC gas. 
Interferometric S-Z measurements, and upcoming high resolution observations 
with bolometric multi-frequency arrays containing large number of elements, 
expand the scope of S-Z science and further motivate more accurate treatment 
of the basic effect and its realistic modeling. 

In this paper we clarify various aspects of a relativistically exact description of 
the S-Z effect, and derive accurate analytic expressions for the intensity change 
and crossover frequency that are valid over wide regions of parameter space. 
Several different approaches have been adopted in the relativistic calculation 
of the effect, and although the approximate analytic expressions that were 
derived in these treatments are either consistent when compared at the same 
level of accuracy, or yield consistent numerical results, it has not yet been 
formally shown that the three main approaches to relativistic calculation of 
the S-Z effect are indeed equivalent. This equivalence is formally proved in 
Section 2.1. In order to provide more accurate analytic approximation to the 
thermal component of the effect, we derive - in section 2.2 - an analytic fit to 
the thermal component of the effect which is accurate to orders r^0*. In view 
of the possibility that in some rich clusters r is as high as 0.03, the approximate 
analytic expansion to first order in r and to eighth order in necessitates also 
the inclusion of multiple scatterings, of order r^. This has been considered by 
Itoh et al. (2001), and is treated further here in section 2.2, where we consider 
the contribution of two successive scatterings (terms to order r^) and find that 
the r^0^ - r^0^ contributions are appreciable near the crossover frequency 
for which we derive an analytic approximation. In section 2.2 we also include 
an analytic expression for the change in the photon occupation number that 
is accurate for gas temperatures as high as 50 keV. In sections 2.3 and 2.4 
the kinematic and the full effect are calculated, including the cross terms that 
depend on both the gas temperature and the cluster velocity. Our results for 
these terms differ somewhat from those obtained by other groups. The last 
section (3) contains a more general discussion of the results presented in this 
paper. 
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2 Exact Description of the S-Z Effect 

In this section we explicitly calculate the thermal, kinematic and cross terms of 
the S-Z effect. While such calculations have already been carried out in various 
papers, our approach here is more direct, and the results are roughly at the 
same level of accuracy as achieved in the recent work of Itoh & Nozawa (2003), 
where the results of a numerical fit are given in a tabulated form. Following 
a detailed description of the full effect, and our analytic approximations, we 
present a fit to the exact calculation that is valid at gas temperature well 
beyond 15 keV. We begin with a general comparison of the main theoretical 
descriptions of the effect. 

2.1 Equivalent Relativistic Treatments 

Three main approaches have been adopted in the calculation of the relativistic 
S-Z effect, those employed by Rephaeli (1995a), Sazonov & Sunyaev (1998), 
and Nozawa et al. (1998a). Different analytic approximations to the exact 
expressions that were obtained in these treatments were found to be consistent 
when compared to the highest common order of ©. Here we formally prove 
that these three approaches are indeed equivalent. Basically, this follows from 
the fact that although the method used by Nozawa ct al. (1998a) is more 
general than those used by Rephaeli (1995a, 1995b) and Sazonov & Sunyaev 
(1998), the simplifying assumptions made in deriving analytic approximations, 
namely neglecting electron recoil, and expansion in powers of a small quantity 
- the ratio of the change in the photon energy to the photon energy - render 
it equivalent to the other two approaches. 

Consider first the method employed by Rephaeli (1995a, 1995b). The probabil- 
ity of scattering of an incoming photon (direction fiQ = cos 6*0) to the direction 
Hq — cos 9'q is (Chandrasekhar 1950) 



where the subscript refers to the electron rest frame. The frequency shift is 



where u, v' are the frequency of the photon before and after the scattering, 
and j3 = v/c is the dimensionless electron velocity in the CMB frame. It is 
somewhat more convenient to use the variables /i, s and (3 instead of /i, /i' and 
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(2) 
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(3. The probability that a scattering results in a frequency shift s is (Wright 
1979) 



e'/ (/ip, ^o ) 



'P{s^P) = ^J :! ''TZ d,,. (3) 



Averaging over a Maxwellian distribution for the electrons yields 



J^Ve-^p;.,.)./. ^ (4) 

The total change in the photon occupation number along a line of sight (los) 
to the cluster, due to scattering off electrons with thermal {t) velocities, can 
now be written as 



Ant{x) J [n {xe') -n{x)] Vi (s) ds (5) 

—oo 

where x is the dimensionless frequency, x = hu/kT, and t = J riedl is the 
optical depth to Thomson scattering. 

In the approach of Sazonov & Sunyaev (1998), which is similar to that of 
Rephaeli (1995a), the starting point is the photon transfer equation 

^(/^o, ^o) = cueo J [no (//o, i^o) - no (//q, i^o)] (6) 
where the Planckian occupation number is 



and Heo & fiQ are the electron number density and the cosine of the angle 
between the directions of the electron motion, and the photon velocity as 
seen in this frame, respectively. The differential scattering cross section in the 
electron rest frame is 



\'^Tf{l^O,IJ''o), (8) 



dQ'o 2 



where ctt is the Thomson cross section. Lorentz transformation of equation 
(6) with Ueo — n^/j leads to 
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where 



5n = no (/Xq, uq) - no (/xq, i^q) = ^ - (10) 



Now, in the Thomson hmit, 



1 + M . . 

and 



dnt f f f da 5n -i^ , i r 02 



dr J J J dn'o (1 + pfXoY 
where 



e-^dfjLodfi'olP^dp/D (12) 



1 



_7-l 



D = J e-^7^/3'd/3. (13) 


This is exactly the expression that was derived by Rephaeh (1995a, 1995b) 
using a different approach, hence the two methods are shown to be fully equiv- 
alent. 

In the third method, that of Nozawa et al. (1998a), the starting point is a 
more basic expression for the time rate of change of the photon occupation 
number given by a generalized Kompaneets equation 



dt^ ^ J J J {2ny 



^{uj)^-2 I I I ^^d^p'd^k'W[n{uj){l + n{uj'))Pe{E) 



-n{u;'){l + n{u;))P,{E')], (14) 



where 



W=l—] S^p + k-p' -k'), (15) 

VP Pj \P pj \P pj 

p = -2 (p • k) = -2uE (1 - Pfio) , (17) 

p' = 2{p-k') = 2uj'E{l-f3f,',), (18) 
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p and k {p' and k') are the ingoing (outgoing) electron and photon 4-momentum 
{p = {E,p), k = {uj,k), p' = {E',p'), k' = {u',k')), and /xq (/^o) is the cosine 
angle between k and p {k' and p') in the electron rest frame, m and e are the 
electron mass and electric charge and Pe{E) is the electron energy distribution. 

In the Thomson hmit, lOq — loq, therefore E — E' , and Equation (14) simphfies 
to 



Now X is essentially the differential cross section, equation (8), used by 
Rcphacli (1995a) and Sazonov & Sunyaev (1998) in a less general form. The 
approximations made by Nozawa et al. (1998a) are based on the fact that 
Ax ^ 1 in the CMB frame (which is in fact, used as a perturbation parame- 
ter), and that T/Te <C 1. However, since hu / {rrieC^Q) = hu/kTe ~ T/Te -C 1, 
it follows that hu/rrigC^ <C O, and the Thomson limit is valid. Note that 
Ax <^ 1 implies that ^ 1, and therefore the condition hv/{mf.(?) <^ 1 still 
holds in the electron rest frame. It can now be seen that even though the treat- 
ment of Nozawa et al. (1998a) is seemingly more general, their final expression 
- in the limit when the commonly adopted approximations are valid - reduces 
to those obtained by Rephaeli (1995a) and Sazonov & Sunyaev (1998). There- 
fore, all the three methods are equivalent and will therefore yield equivalent 
results to all orders in 0, provided the expansion in is valid (0 <C 1). 



2.2 The Thermal Component 



2.2.1 Single Scattering Limit 

In an optically very thin medium with r < 0.01, only single scatterings of 
CMB photons need to be (practically) considered, and the leading term in the 
expression for An is (by far) of order r. Using equation (9) with an integral 
over (relativistically generalized expression for) a Maxwellian distribution of 
electron velocities, we obtain 



111 

dt D J J i ^ 7^ (1 + PfJ^o) d% / M M A-o w y J 



(3=0 ^^=-l//o= 

To first order in r. An is 



7 



To perform the integrations in this expression, wc expand the integrand in 
powers of /3 around (3 — with the derivatives (at P — 0) compactly written 
by defining 



A„ = -2"-2x"[sinh-' ^]("-^), (22) 

where /*^"^ denotes the n'th derivative of /. Carrying out the expansion to 
and integrating yield 



— ^ = OFi + e'Fs + e^Fs + e^F4 + e^Fg + e^Fg + e^Fy + e'Fg (23) 

T 

where Fj are listed in Appendix A. 

It is instructive to compare our result in equation (23) to equations (2.25)- 
(2.30) of Nozawa ct al. (1998a) that were obtained using a different approach. 
The functions Ai can be expressed in terms of X and S 



X = xcoth — , 

and the relations between A^, X and S are given in Appendix B. Substituting 
the expressions in Appendix B back into Fi - F5 and equation (24), we obtain 
exactly equations (2.25)-(2.30) of Nozawa et al. (1998a) term by term. This is, 
of course, just what we expect given the proven equivalence of the approaches 
discussed above. 



2.2.2 Multiple Scatterings 

Although formally very accurate, the single scattering analytic approximation 
to order 0^ is actually limited to small values of r. This optical depth is 
typically ~ 0.01 in a rich cluster, but values as high as of up to 0.03 are not 
unexpected in some clusters. Thus, the possibility of double scatterings is not 
ignorable to the level of accuracy afforded by the above expansion to order 
r©^. A more self consistent calculation must include the higher order terms 
up to r^G^, which can be comparable to the r0^ term. The importance of 
considering the case of finite optical depth was already pointed out by Molnar 
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& Birkinshaw (1999), Itoh et al. (2001), and Dolgov et al. (2001). Our goal 
here is to extend the analytic approximation by calculating all the relevant 
contributions up to and including terms that are of order r^G^. 

The impact of second scatterings on the occupation number can be calculated 
iteratively by describing the incident radiation in terms of the occupation num- 
ber, 111, that results after the radiation has already been singly scattered. In 
this way one can extend the treatment to the case of any number of scatterings 
(Wright 1979). Now, rii is simply 



ni = no + Ano, (25) 

where Ano is given in equation (23). Substituting ni into the right hand side of 
equation (12) we can solve for Ani. Note that the result of the naive iteration 
should be divided by 2! to avoid double counting; a given total frequency shift 
{ly'/v — z) results from two scatterings irrespective of the order of the two 
scattering events. The expression for the occupation number of the doubly 
scattered radiation is 



n' = (no + Tno®V+ ^t^uq ®V , (26) 

where ® denotes convolution and (as before) V is the scattering probability, 
equation (3). Expanding the exponent to second order in r, we have 



Ant = - ® ^ ~ ^) + 2"^^ *^ ® ~ ^ ® ■ 

The first term on the right hand side of equation (27) {n^V—n) was evaluated 
in section 2.2.1; here we calculate the contribution. Performing the integra- 
tions in a way similar to the calculation of the first order r contribution, we 
obtain 



+ T'e^Fi2. (28) 

Altogether, the resulting change in the occupation number calculated to orders 
tG^ and r^G^ is 



An* = t(GFi + G'Fa + G^Fs + G^F4 + G^Fs + G^Fg + G^Fy 

+ G^Fg) + (G^Fg + G^Fio + G^Fn + G^Fis) . (29) 

It can be easily shown that equation (29) satisfies photon number conservation. 
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The terms are important especially near the crossover frequency, whose 
vahie depends both on and r. as derived expUcitly below. Also of interest 
is the unignorable effect of double scatterings at frequencies on the far Wien 
side of the spectrum, due to the fact that the total change in the frequency 
of the photon after two scatterings is particularly high. More specifically, we 
examined the contribution of these terms to A/. For r in the range [0.01,0.04] 
and in the range [1/130, 1/35], the corrections due to these terms are larger 
than 1% near the crossover frequency as expected. Clearly, the spectral band 
around this frequency - for which the corrections are not ignorable - widens 
as increases. 



2.2.3 Analytic Fit to the Thermal Effect 

The observed range of IC gas temperatures is ~ 3 — 15 keV, when the gas is 
taken to be isothermal (in the fit to the observed X-ray spectrum). High spatial 
and spectral resolution measurements with the Chandra and XMM satellites 
show more and more that the temperature distribution is non-uniform, which 
is only to be expected given the hierarchical growth of structure in the con- 
text of CDM models. A more realistic modeling of IC gas yields hotter and 
colder regions in some clusters, perhaps reficcting the nature of the forma- 
tion process of clusters from aggregates or sub-clumps. In a naive modeling 
of the temperature evolution of clusters a virial mass-temperature relation is 
sometimes obtained with clusters of a given mass having a higher temperature 
the earlier they virialized. If so, a roughly linear relation would be predicted 
between the temperature and the virialization rcdshift of the cluster. Irrespec- 
tive of the reality of such a scaling, it is of interest to extend the accuracy of 
the analytic approximations for An to values of kT^ which are higher than 15 
keV. Equation (23) is not sufficiently accurate at higher temperatures because 
the series expansion docs not converge and rapidly oscillates around the exact 
numerical result (equation 21). 

Here we obtain a more accurate fit to the results from the exact calculation 
using a polynomial function of the form 

5 in 5 

An(x)/r = ^ ^ 6,, A,0' = ^ (30) 

i=l j=l i=l 

where the numerical coefficients bij are obtained by means of a best fit to 
the exact expression (equation 21). The coefficients are given in Appendix 
A. This expression is in agreement with the exact calculation to better than 
2.5% for < 1/10, or kT^ < 50 kcV. This level of accuracy is higher than 
that obtained by Itoh et al. (2001) who used the analytic expansion (equation 
23). Evidently, higher accuracy is attained by not insisting on using the exact 
form of the analytic expansion also at higher temperatures; our approach is 
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thus functionally less constrained. Very recently, Itoh & Nozawa (2003) have 
extended their previous analytic approximation for An to higher temperatures 
by fitting the results of the exact numerical calculation with a polynomial of 
order for a; < 1.2, and a fully numerical fit to a polynomial in powers of 
both © (to 13th order) and x ior x > 1.2 (with the large array of coefficients 
given in a tabulated form). This more cumbersome procedure does result in 
a more accurate representation of the exact numerical calculation, which is 
generally within 0.4%. 

2.2.4 Crossover Frequency 

As was pointed out in Rephaeli (1995b), Sazonov & Sunyaev (1998) and 
Nozawa et al. (1998a), the crossover frequency (xq) depends on the gas tem- 
perature. In previously derived expressions for the functional form of Xq, only 
the dependence on © was considered, and deduced from a fit. In order to 
find an analytic approximation to the shifted crossover frequency we write 



and substitute this expression in equation (29). Expanding in the small param- 
eters T and ©, we then obtain the values of the coefficients in the expression 
for xq by setting Arif — 0, 



It can be readily shown that, as expected, Xq monotonically increases with r 
and ©. 



2.3 The Kinematic Component 

In an idealized cold IC gas all electrons have the same velocity in the CMB 
frame - the cluster peculiar velocity, Pc- Substituting ioi (3 in equation (9), 
and applying the Lorentz transformation 



Xq = 3.8300161 (1 + Ci© + C2©^ + C3©^ + 04©"^ + C5T + CqtQ 




(31) 



Xq = 3.830016 (1 + 1.120594© + 2.078258©^ - 80.748072©=^ 
+ 1548.250996©^ + 0.800424r© + 1.183073t©^) . 



jJ'c - Pc 

1 - /5cAic 



(32) 
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where fic is the angle cosine with respect to the los of an observer (as measured 
in the CMB frame), we obtain 



^n=—T{l-l3clIc) 



></(3-,?) + K->)(f^) 




(33) 



where 5n is given in eq. (10), and 



z — 



1 + Pf^o 

1 + Pflo 



7e'(l+«)(l-/3c/^c). 



(34) 



We evaluate the latter expression for An as part of the calculation of the full 
effect in the next section. 



To first approximation, when and /3c are sufficiently small, the impact on 
the CMB by Compton scattering in a cluster with finite peculiar velocity is 
simply a sum of the thermal and kinematic components. For higher values of O 
and j3c, higher order terms have to be included. As a result, cross terms linking 
the two components arise. The inclusion of the peculiar velocity of the cluster 
in a consistent calculation (not as a separate effect) was originally considered 
by Sazonov & Sunyaev (1998), and later also by Nozawa et al. (1998b). Their 
calculations yield consistent results for the cross terms that depend on both 
and f3c- We repeat the calculation of the full effect by generalizing our 
treatment described in the previous sections. We do so in part because our 
calculation yields slightly different results for the cross terms, although our 
results for the pure thermal and kinematic components are exactly identical 
to those obtained in the latter papers. 

We begin with the transformation of the electron thermal energy and momen- 
tum, Et and pt, from the cluster to the CMB frame, E and p, respectively 



Applying these transformations in the expression for the Boltzmann factor we 
obtain 



2.4 The Full Effect 



E = -fc{Et + c(3c-pt) 




(35) 
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exp 



exp 



7t-l 

e 



exp 



7c7(l-/9-/9c)- l 

e 



(36) 



where (3 and 7 are the velocity and Lorentz factor in the CMB frame. 
Since (Pp/E is invariant under Lorentz transformations, 



E 

d^Pt = -^d?P = 7c (1 - /3c ■ 0) p^dpdi^dip 
^^,[l-P,-P)^'(3''d(3diid^. 

Expanding to second order in (3c we obtain 



(37) 



7t - 1 



p^dpt exp ^ 
= 7c (1 - /3c ■ pt) exp 



7c7(l-/3-/3c)-ll„,,, 



e 



X exp 
Note that 



l + /3c-/3(|-l^l + 



e 



2 
e 



(38) 



(3c- (3 — (3c(3 (cos 6 cos ^c + cos (/? sin ^ sin ^c) 

cos 6q + (3 



cos9 — 



1 + (3 cos ^0 



(39) 



where 9 and 9c are angles in the cosmic frame and 9q is in the electron rest 
frame; we select axes such that v?c = 0). Integrating over the azimuthal angle 
(p (divided by 27r) and expanding again near /3c = 0, we obtain 



pidpt exp 



It -I 

e 



1 + 



Afc.(|-l) + ift=/J=(.^-i)(3."-l) 



2^,2 



662 



30 2 20 



exp 



7-1 

e 



7^/3^ci/3. 



(40) 
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Substit.Titing this into the transfer equation (Equation (21)) and performing 
the integrations yield the correction to the thermal effect, Equation (23) 

+ p^^ (^F, + eF22 + e'F23 + e=^F24 + e^Fss) , (4i) 

where Fj are found in Appendix A. 

The purely kinematic terms in equation (41) are identical to the corresponding 
terms in similar equations obtained by Sazonov & Sunyaev (1998) and Nozawa 
et al. (1998b). However, our results for the cross terms differ somewhat from 
those obtained in the latter two papers. The difference is due to integration 
of the transfer equation over '~f^f3'^df3dfj,d(f. rather than over equation (37) - 
as we have done here - which includes the extra factor 7c(l — /5c • /5). From 
a practical point of view this difference is quite small, typically less than 1%. 
We note that both our result, equation (41), and the results derived in all the 
latter three papers satisfy photon number conservation. 



3 Discussion 

In this paper we have first shown that the three main, seemingly different 
approaches to the calculation of relativistic corrections to the expressions of 
Sunyaev & Zeldovich (1972), are indeed equivalent in the limiting case for 
which they are derived. We then extended the accuracy of the analytic ap- 
proximation to the thermal and kinematic terms such that the purely thermal 
component of the effect is accurate to terms of order B^, and the cross terms 
to order /^c©^- We have also obtained a relatively compact fit to the exact 
numerical calculation of An that is accurate to within 2.5% for gas temper- 
atures as high as 50 keV. Either this functional fit, or the more accurate - 
but perhaps somewhat less convenient - fully numerical tabulation of Itoh & 
Nozawa (2003), can be used in the analysis of SZ measurements - including 
cosmological parameter determination - whose current observational errors 
are well above the accuracies provided by these fits. 

While the accuracy of these expressions is limited very close to the crossover 
frequency, this is usually of little practical consequence because of the typi- 
cally broad detector spectral response: Even when the measurements are made 
very close to the crossover frequency for the purpose of measuring cluster pe- 
culiar velocities, what is required is only a sufficiently accurate value of the 
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AI convolved over the spectral response. The analytic expression given here 
is accurate to within a few %. Convolving both the exact (eq. 21) and the 
approximate expressions (eqs. 23 & 30) over a Gaussian detector response 
centered at 214 GHz, and with a FWHM of 30 GHz (characteristics of the 
second MITO channel), we find that the two results agree to within 0.2%. 
Thus, the impact of this small difference (close to the crossover frequency) on 
the determination of cluster peculiar velocities is expected to be negligible. 

The accuracy of the theoretical description of the SZ effect surpasses the pre- 
cision with which the effect can be measured at present. In addition to ob- 
servational errors, systematic modeling uncertainties are relatively large. One 
of these uncertainties stems from the unknown gas temperature profile. From 
a theoretical point of view this can be readily assessed for a given assumed 
profile. For example, consider the case of an idealized (and perhaps somewhat 
unrealistic) polytropic equation of state P oc p^, where 7 is a free parameter 
which is determined from X-ray measurements. In order to generalize our re- 
sult for the analytic approximation for An to the case of non-isothermal gas, 
we need to replace the spatially dependent parameters by their appropriately 
weighted expressions. Consider, for example, the commonly used /3 model for 
the gas density, n^ir) oc [1 -|- (r/rc)^]"^'^'^^, where Tc is the gas core radius. 
(Note that here 7 and j3 denote different quantities than those in the rest of 
the paper.) Since Tg oc p^^'^, the products of r0* that appear in the analytic 
approximation for An have to be replaced by weighted averages along the los 



where z is the radial extent of the gas distribution in units of the core radius, 
no & To are the central values of the electron density and temperature, and 
To = 2aTnorc. 



Appendix A 

The following is a list of the functions Fj, Gi and Hi in terms of the functions 
Aj [defined in Equation (22)] 
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(42) 



—2 



Fi = 2A, + ^A2 

47 21 

F2 = 5Ai + —A2 + —A3 + 




11 



1920 



15 



+ 



15 , 2505 , 3549 , 14253 , 9297 , 

F4 = Ai H A2 H A3 H H A5 

^ 4 32 ^ 20 ^ 160 560 ^ 

12059 1 1 

+ "8960"^^ + 21^^ + 1680^^ 
^ 135 , 30375 , 62391 , 614727 , 124389 , 

355703 2071 1879 11 11 

^120"^^ + + 6720^^ + 1792'^'' + 215040^^° 

45 , 7515 , 28917 , 795429 , 2319993 , 

— Ai A2 H A H Aa H A^ 

16 128 16 128 448 

12667283 , 201631 , 10655 , 46679 , 10853 , 

Aq H A7 H As H Ag H Aio 

7168 672 384 32256 258048 

29 , 29 , 

H All H ^12 

46080 7741440 

7425 128655 360675 50853555 45719721 

- "^^^ + ^09^^^ + ^5^^^ + 2048 + 1024 
458203107 22251961 71548297 26865067 

+ 16384 + 2688 + 53760 + 215040 

7313155 , 1123 , 6361 , 37 , 37 , 

H Aio H All H A12 H H Au 

1032192 4608 1290240 691200 154828800 

^ 675 , 6345 , 86751 , 28579473 , 463090581 , 

^« = 16^^ + 12^^^ - + -^^^ + ^79^^^ 

8680356807 , 407333911 , 304758409 , 14281971623 , 

H A H Ay H Ag H 

28672 2688 7680 2365440 

32154229291 , 17138321 , 36841447 , 481771 , 

H ^10 H ^11 H ^12 H ^13 

56770560 506880 28385280 15206400 

229693 ^^4 

^ 486604800 ^ 5913600 ^ 1703116800 

^ . . 17 . 3 , 1 , 

F9 = 4Ai + — ^2 + + ^^4 

, 295 , 873 , 1271 , 21 7 , 

Fio = 20yli H A2 H A3 H A4 H A5 H Ag 

10 1 ^ ^ ^ 20 ^ 160 40 640 

44643 , 19557 , 27549 , 34873 , 

Fii = 40Ai + 599^2 H A3 H A4 H A5 H Aq 

50 50 400 6400 

367 , 367 , 

H A7 H As 

1920 153600 

45 , 96651 , 8659449 , 62384943 , 38586081 , 

F12 = — Ai H A2 H A3 H A4 H A5 

2 32 800 6400 11200 

103117227 , 82813 , 590831 , 859 , 859 , 

H Ae H Ay H As H Ag H Am 

179200 1680 268800 17920 2150400 

9 47 7 

^13 = T^A + — A2 + —A3 



Z 20 4U 

35 , 1023 , 833 , 183 , 11 , 

Fi4 = — Ai H A2 H A3 H A4 H A5 

4 ^ 40 80 ^ 160 320 ^ 

175 , 18721 , 45969 , 53757 , 

Fi5 = Ai H A2 H A3 H A4 + 

16 160 320 1120 



53757 13291 3917 
1120"^' + ^240"^' + 13440^' 
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210 ^ 

315 , 45231 , 324387 , 355815 , 1020813 , 
Fia = -^^1 + + -^As + ^A, + -^A, 

781943 1983 167 11 
+ ^6880"^' + 1120^' + 3360^^' + 21504^' 
1 . 11 . 
^"=3^' +120-*^ 

4 , 169 , 3 , 19 , 

"3 60 4 480 

10 , 2609 , 2349 , 31277 , 1227 , 23 , 
Fi9 = — Ai H A2 H A3 H A4 H A5 H A 

3 120 120 6720 3360 2688 

5 , 14957 , 34113 , 223653 , 14963 , 16249 , 
F20 = -Ai H A2 H A3 H A4 H Ak H Aq 

" 2 160 160 1792 560 6720 

941 ^ 101 , 

H A7 H As 

10080 80640 

205 , 81869 , 6291 , 6482541 , 2372827 , 

F21 = Ai H A2 H A3 A4 A5 

64 512 160 8960 4480 

6820873 , 269683 , 109183 , 515 , 1 , 

Ae Ay Ag Aq Aio 

53760 20160 161280 32256 7168 

8 . 23 , 7 , 7 , 

F22 = -Ai H A2 H A3 H A4 

22 31^6 ^ 10 ^ 240 

20 , 367 , 427 , 1967 , 11 , 11 , 

F23 = — Ai H A2 H A3 H A4 H A5 H Aq 

^ 3 12 20 480 40 1920 

2131 , 19369 , 76163 , 3099 , 12059 , 

F24 = 5Ai H A2 H A3 H A4 H A5 H Ae 

24 1^ 16 ^^ 80 640 140 6720 

4 . 1 . 
H A7 H Ag 

63 1260 

4925 , 126903 , 1085261 , 185183 , 

F25 = -5Ai H A2 H A3 H A4 H A5 

^ 16 80 640 280 

522981 10387 9403 55 11 

+ ^480"^' + 1008"^' + 20160^' + 5376^' + 129024^^° 

Hi = 2Ai + 7A2 
4 

H2 = 4.6699698Ai + 5.4022055A2 + 0.77091691A3 - 0.040304645A4 

- 0.0138114993A5 - 0.0012957652A6 - 7.168369 x lO'^Ay 

- 2.3560827 x 10"% - 4.333988 x 10"% - 3.4497085 x lO^^^Aio 

H3 = 19.0998386Ai + 57.853156A2 + 39.0362A3 + 9.6917067A4 

+ I.2O6357A5 + 0.09119917A6 + 0.004517831A7 + 1.4162735 x 10"% 
+ 2.5544641 x lO'Mg + 2.0823368 x lO^^Aio 
= 349.87908Ai + 249.46964A2 + 105.99813A3 + 41.012553A4 
+ 8.7994039A5 + 0.96698562A6 + 0.058225778A7 + 0.0019859754A8 

Uo 

H5 = -2818.9102Ai - 3IO9.3757A2 - 1577.9854A3 - 465.38533A4 

- 79.OO9596A5 - 7.7208374A6 - 0.44258019A7 - 0.01482893A8 



+ 3.59624 X 10"^ Ag + 2.6060166 x 10"% 
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- 2.6904478 x 10~% - 2.0134634 x 10~^Aio 

(A.l) 



Appendix B 

From Equations (22) and (24) the following relations are obtained 



^.sinh(|)=-i5 

A2 sinh = SX 

A4 sinh = ASX^ + SS^X 

A5 sinh = -8SX^ - AAS^X^ - 8S^ 

Ae sinh (^0 = 16^1^ + 208S^X^ + 136S^X 

Aysinh (^1^ = -32SX^ - 9125^X^ - 1440^^^^ - 1365^ 

X 

2 

X 

2 



As sinh f I j = 64^X^ + 3840^^X^ + 12288^^1^ + 3968^^X 
Ao sinh ( ^] = -128SX^ - 158Q8S^X'^ - 929285^1^ - 686085^X2 



- 3968^^ 
= 256^X'^ 

+ 176896,S''X. 



^10 sinh f = 256^X'^ + 64256^^1^ + 652416^^^^ + 920576^^^^ 



(B.l) 
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Deviations of Approximations From the Exact Calculation (15keV) 
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Fig. 1. The deviations of the fifth (Itoh et al. 2000a) and eighth order approximations 
(equation 23) from the exact result. 
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